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We develop the theory of light propagation under the conditions of electromagnetically induced 
transparency (EIT) in systems involving strongly interacting Rydberg states. Taking into account 
the quantum nature and the spatial propagation of light, we analyze interactions involving few- 
photon pulses. We demonstrate that this system can be used for the generation of nonclassical 
states of light including trains of single photons with an avoided volume between them, for imple- 
menting photon-photon quantum gates, as well as for studying many-body phenomena with strongly 
correlated photons. 
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The phenomenon of electromagnetically induced trans- 
parency (EIT) [1 in systems involving Rydberg states 
[2] has recently attracted significant experimental [3HIU] 
and theoretical [TTH2T] attention. While EIT allows for 
strong atom-light interactions without absorption, Ryd- 
berg states provide strong long-range atom-atom inter- 
actions. Therefore, the resulting combination of EIT 
with Rydberg atoms is ideal for implementing mesoscopic 
quantum gates [2j [16] and for inducing strong photon- 
photon interactions, with applications to photonic quan- 
tum information processing [2j [TTHT4| [T9H22] and to the 
realization of many-body phenomena with strongly inter- 
acting photons [23]. At the same time, the many-body 
theoretical description of EIT with arbitrarily strongly 
interacting Rydberg atoms, taking into account the full 
quantum dynamics and the spatial propagation of light, 
has not been reported previously. 

In this Letter, we develop such a theory by analyzing 
the problem for at most two incident photons, which, in 
turn, provides intuition for understanding the full multi- 
photon problem. We show that Rydberg atom inter- 
actions give rise to photon-photon interactions, which, 
below a critical inter-photonic distance, turn the EIT 
medium into an effective two-level medium. This can 
be used to implement photon-atom and photon-photon 
phase gates and to enable deterministic single-photon 
sources. 

The basic physics is illustrated by considering a simple 
case [Fig. gb)], in which a single-photon wavepacket £ 
propagates in an EIT medium [level scheme in Fig. [TJa)] 
with a central control atom at z = prepared in a Ry- 
dberg state \r'). Atoms in another Rydberg state |r), 
coupled by the EIT control laser [Fig. ma)], experience 
a van der Waals potential V(z) = Cq/zt due to the in- 
teraction with the control atom, which is decoupled from 
the applied fields. Alternatively, one could apply an elec- 
tric field to induce dipole moments in states |r) and |r'), 



resulting in V oc 1/z 3 [2] . 

Sufficiently far away from z = 0, the incident photon 
propagates in a standard EIT medium featuring a two- 
photon-resonant control field with Rabi frequency Cl. In 
the vicinity of the control atom, however, the state |r) 
is shifted so strongly out of resonance that the photon 
sees only a two-level (|g), |e)) medium with transition 
linewidth 27. The critical z, at which the interaction is 
equal to the EIT linewidth, separates these two regimes 
and corresponds to the Rydberg blockade radius [TTJ [24] . 
When the single-photon detuning A = 0, the resonant 
blockade radius z\> is thus defined by V(z\>) = Ct 2 /j 
(h = 1), while for A ^> 7, we define the off-resonant 
blockade radius zb via V(zb) = ^ 2 /A (we assumed 
A/Cq > 0). Since the blockade region extends over 2z b ( B ) 
[Fig. [IJb)], the presence of the control atom locally cre- 
ates an absorbing or refractive medium with an effective 
optical depth db(B) = 2cbb(B)/£? where d is the resonant 
optical depth of the \g) — \e) medium (Cl = 0) of length L. 
Interesting effects occur at large blockaded optical depths 
c?b(B)- I n the resonant case, assuming ^> 1, the pres- 
ence of the \r') excitation causes complete scattering, i.e. 
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FIG. 1. (a) EIT level scheme, in which a ground state \g), an 
excited state |e), and a strongly interacting Rydberg state |r) 
are coupled by a quantum probe field S and a classical con- 
trol field with Rabi frequency Cl and single-photon detuning 
A. (b) Interaction of one photon with a Rydberg excitation 
stored at z = 0, which modifies the propagation within the 
blockade region \z\ < Zb(B)- (c,d) Interaction of two counter- 
propagating (c) or co-propagating (d) photons. 
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loss, of the incoming photon. Off resonance, for ^> 1 
and 6^(7/ A) 2 <C 1, the interaction with the Ir^-atom 
imprints a phase ~ g?b7/A on the probe photon and re- 
duces its group delay by ~ ^b7 ? as its group velocity 
is increased to the speed of light, c, within the blockade 
region. 

In the off-resonant case, this simple system has di- 
rect practical applications. First, by encoding a qubit 
in the ground and \r') states of the central control atom, 
one can implement a phase gate between the probe pho- 
ton and the atom. Second, the protocol of Ref. [25] al- 
lows to implement a phase gate between two photons 
by successively sending them past the control atom that 
is appropriately prepared and manipulated between the 
passes. Selective manipulation of the control atom can be 
achieved particularly simply if it is of a different species 
or isotope. Third, a phase gate between two photons can 
also be achieved by storing one of them in the \r') state 
of the control atom and sending the other one through 
the medium. While storing a single photon in a single 
atom is difficult, the same effect can be achieved by stor- 
ing [26] [27] the photon in a collective \r') excitation, as 
we will discuss below. 

The results of this simple problem can be extended to 
the case of multi-photon EIT propagation in Rydberg 
media. First, off-resonance, two counter-propagating 
photons [Fig. [ljc)] can pick up a phase ~ g£b7/A, en- 
abling the implementation of a two-photon phase gate 
[T2j [14] . Second, a pulse of co-propagating photons [Fig. 
[TJd)] will evolve into a non-classical state corresponding 
to a train of single photons [19] and exhibiting correla- 
tions similar to those of hard-sphere particles with ra- 
dius Zb(B)/2. These correlations arise from scattering of 
photon pairs within the blockade region. Third, in the 
regime where z^ is larger than the EIT-compressed pulse 
length, a, both co- and counter-propagating resonant se- 
tups might be usable as single-photon sources since all 
but one excitation will be extinguished. In the follow- 
ing, we present a detailed theoretical analysis of these 
phenomena. 

Interaction of a photon with a stationary excitation. — 
We begin by detailing the solution of the problem of a sin- 
gle photon propagating in a medium where state |r) expe- 
riences a potential V(z) [Fig.JTJb)]. Treating the medium 
in a one-dimensional continuum approximation, working 
in the dipole and rotating- wave approximations, and adi- 
abatically eliminating the polarization on the \g) — \e) 
transition, the slowly varying electric field amplitude £ 
of the single-photon wavepacket and the polarization S 
on the \g) — \r) transition obey [26] [27] 



(d t + cd z )£(z,t) 

d t S(z,t) = -iU(z,t) - l -fS(z,t) 
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a^S(z,t), (1) 



£(z,t). (2) 



Here T = 7 — zA, U(z,t) = V(z)S(z,t), g is the atom- 
field coupling constant, and n is the atomic density. We 



have neglected the depletion of state \g) and the finite 
lifetime of the Rydberg state |r), which is typically much 
longer than the propagation times considered here [2]. 
Assuming that all atoms are in state \g) before the arrival 
of the photon, Eqs. (T|2) can be solved to give 



where 



duje 



dzx(z,i 



kx(z,oj) 



d~f[uj - V(z)] 



Lfi 2 - (A + «7)[w- V(z)} 



(4) 



and £(— L/2,cj) is the Fourier transform of the 
wavepacket incident at z = —L/2. 

For narrowband pulses, we expand x 111 00 an d, assum- 
ing A ^> 7 and L ^> 2^b, reduce Eq. ([3| to 



£ (L/2, t)&£ (-L/2, t - L'/v e ) e^, 



(5) 



where v g « cQ 2 /(g 2 n) = 2Q 2 L/(d>y) is the EIT group 
veclocity. In order to avoid the Raman resonance at 
V + Q 2 /A = 0, we assumed A/Cq > 0. Since the photon 
travels at c within the blockade region, the group delay 
comes from a reduced medium length L' = L — ^ttzb ~ 
L — 2zb- Additionally, the emergence of a two-level 
medium within \z\ < zb gives an intensity attenuation of 
e~ 2r > with 2rj = ffd B (7/A) 2 « ^(7/A) 2 and a picked- 
up phase of ip = — ^^(7/ A) ~ — ^g?b(7/A). Thus, with 
g£b ^> 1 and a properly chosen A ^> 7, one can get a 
considerable phase and/or change in group delay with- 
out significant absorption. For the same derivation on 
resonance (A = 0), the main effect is an intensity atten- 
uation of « exp(— db)j as expected for a two-level medium 
of length 2^b- 

It is straightforward to extend our analysis to a delo- 
calized \r') excitation, i.e. a spin wave, that is spread over 
many atoms. Far off resonance, the effect of the control 
atom is independent of its position, such that a single 
control atom and a corresponding spin wave affect the 
incident photon identically. On resonance, with d^ ^> 1, 
the \r') spin wave causes complete scattering of the in- 
coming photon. At the same time, after tracing out the 
scattered photon, which carries information about the lo- 
cation of the scattering, the spin wave itself disentangles 
into a classical mixture of pieces of length ~ z\>. 

Interaction of propagating photons. — We now consider 
the problem of propagating photons interacting with each 
other. Regarding £ and S in Eqs. ( |1|2[ ) as operators 
with same-time commutation relations [£(z), £^{z')] = 
[S{z),S\z f )\ = 5(z - z') [27 and taking U(z) = 
J dz'V(z-z , )S\z')S(z , )S(z) 1 Eqs. {IJ2} become Heisen- 
berg operator equations [28] for the case of photons co- 
propagating in a Rydberg EIT medium [Fig. [TJd)]. Al- 
ternatively, for the case of two counter-propagating pho- 
tons [Fig. [ljc)], we define operators £1(2) an d ^1(2) f° r 
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the right- (left-) moving photon. For Si, the interaction 
is U(z) = J dz f V{z-z f )S\{z f )S 2 {z f )S 1 {z), and vice versa 
for S 2 . 

Since the physics of two counter-propagating pho- 
tons is similar to the spin- wave problem above, 
we begin our analysis with this case [Fig. [TJc)] . 
Letting \^(t)) be the two-excitation wavefunction 
[29], we define ee(z 1 ,z 2 ,t) = <0|^i (^i)^ 2 (^2) |^(^)) , 
es(z 1 ,z 2 ,i) = (0\£ 1 (z 1 )S 2 (z 2 )\^(t)), se(z 1 ,z 2 ,i) = 
(0\Si(zi)£ 2 (z 2 )\i/j(t)), and ss(z u z 2 ,t) = 
(0|Si(zi)S 2 (£ 2 )|^W). Eqs. ([Tj2| then yield a sys- 
tem of equations for these four variables. Defining 
es± = (es ± se)/2, one finds that es_ is small and 
does not significantly affect the dynamics. Dropping 
es_, defining center-of-mass and relative coordinates 
R = {z\ + z 2 )/2 and r = Z\ — z 2 , and taking a Fourier 
transform in time, one obtains c<9 r v = M(r, cj)v, where 
v = {ee(R, r, cj), es+(R, r, u)} and 
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R enters only through boundary conditions and is, thus, 
not important in the present case. For narrowband 
pulses, we can expand M(r, uo) ~ Mo(r) +wMi(r), with 
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Here we defined the effective potential V = TV/ (TV — 
i2Vt 2 ). Outside (inside) the blockade region, V « 
iTV/(2Q 2 ) (V « 1). For \r\ > z b(B) , the two pho- 
tons propagate as dark-state polaritons [26], i.e. we have 
es+/ee = —g^/n/Vt, which is an eigenstate of Mo with 
eigenvalue 0. Since g^Jn ^> £7, the group velocity can be 
read out from the last entry of Mi, which gives twice the 
EIT group velocity v g since the two polaritons propagate 
towards each other. Within the blockade radius, where 
V ~ 1 and V/V ~ 0, the polariton solution ceases to be 



an eigenstate of M , and Eq. (7b) predicts a speed up to 
~ c. Since the time ~ ^b(B)/ c it takes to cross the block- 
ade region is much less than the inverse width of the EIT 
window, the dynamics is highly non-adiabatic such that 
the main result of the interactions is a picked-up factor 
of exp [— J drg 2 nV(r)/(cT)] = ex.p(i(p — rf). This is a 
generalization of the result of Refs. [12j [14] beyond the 
perturbative regime. 

On resonance, 2r] ~ d^. Thus, analogously to the spin- 
wave problem above, the entire EIT-compressed two- 
particle wavefunction decays provided it fits inside the 
medium and d\> ^> 1. The resulting state is a statistical 
mixture of right- and left-moving excitations. 



Off resonance, es + picks up ip \ 
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FIG. 2. (a)-(d) Two-photon counter-propagation for A = 
2O7, Q = 2A, g^fn = 20A and z B = 0.055a, where a is the 
compressed pulse length inside the medium. The color coding 
shows the local phase of ee, while the opacity reflects the two- 
photon density |ee| 2 . The dashed lines are \zi—z 2 \ — zb- The 
full movie is provided in the supplementary material [30] . (e) 
Numerically obtained phase shift ip [we plot cos <p] and atten- 
uation e _r? as a function of cfe compared to the analytical pre- 
dictions (lines). The numerical data corresponds to two dif- 
ferent parameter scans g^/n = 400 A, zb = 0. 0025a, 0.03a 
(dots) and zb = 0.03a, g\fn — 80, 390 (squares). 



diagonal terms in Mo result in a small admixture of the 
bright-state polariton [26], which decays after the wave- 
function exits the blockade region. 

To verify these conclusions, we show in Fig. [2] and in 
the supplementary movie [30] the results of numerical so- 
lutions of the full equations for ee, es, se, and ee in the 
off-resonant case. Despite the bright-polar iton-induced 
oscillations of ee inside and near the blockade region [30 , 
the final phase of the outgoing two-photon pulse perfectly 
agrees with our analytical prediction [Fig. J2je)]. While 
also showing good agreement with the analytical result, 
the obtained loss is slightly larger due to the bright- 
state polariton admixture, which was neglected within 
the above approximate treatment. 

Provided the EIT-compressed two-particle wavefunc- 
tion fits inside the medium, this process, thus, allows for 
the implementation of a nearly lossless phase gate be- 
tween two photons. Taking a specific example of cold 
Rb atoms with |e) = 5 2 Pi/ 2 and \r) = 70 2 Si/ 2 and us- 
ing £1/2tt = 2MHz and A = 2O7, we find zb = 15/im, 
which, for a dense cloud with n = 10 12 cm -3 , gives 
d& = ^X 2 (2zb)ti ~ 9. This yields a significant phase 
of cp « —0.2 and a very small attenuation 2r] ~ 0.02. 
One can increase d^ further by using photonic waveg- 
uides [3TU33] and working with a BEC [34]. 

In the co-propagating case, we define ee{z\,z 2 ,€) = 
(0|f(zi)f(^)|^W), es( Zl ,z 2 ,t) = (0\S(z 1 )S(z 2 )m)), 
and ss(z u z 2 ,t) = (0\S(z 1 )S(z 2 )m)) [Fig. Qd)]. Defin- 
ing es±(z 1 ,z 2 ) = [es(z 1 ,z 2 )±es(z 2 ,zi)]/2, dropping es_, 
and taking the Fourier transform in time, we obtain 
c8rv = 2M(r, u))v. That is, the only difference from 
the counter-propagating case is the replacement of d r 
with (1/2)8r. The resulting equations can be solved 
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A/Cq < 0, the effective potential shows a resonant fea- 
ture, which can give rise to two-polariton bound states. 

We thank S. Hofferberth, M. Bajcsy, T. Peyronel, M. 
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ported by NSF, the Lee A. DuBridge Fellowship, DFG 
through the GRK 792, CUA, DARPA, and the Packard 
Foundation. 



FIG. 3. Time evolution of | ee | 2 for two co-propagating pho- 
tons for Q = 7, g^/n = IOO7, and Zh = 0.08a. The dashed 
lines are \z± — ^2 1 — 2t>, in agreement with the numerical re- 
sults, which show the decay of ee within the dashed lines. The 
full movie is provided in the supplementary material [30] . 



separately at each r. As before, outside the block- 
ade radius, the two-photon dark-state polariton propa- 
gates with group velocity v g . Inside the blockade ra- 
dius, Mq results in fast attenuation on a lengthscale 
~ aK^ 2 ^ 2 )h 2 - This is confirmed by our numerical 
calculations, shown in Fig. [3] and in the supplementary 
movie [30]. Therefore, the two-excitation wavefunction 
evolves into a statistical mixture of a single excitation 
and a correlated train of two photons separated by Zb(B) • 
On resonance, if the photon is scattered when the EIT- 
compressed pulse of length less than is fully inside the 
medium, we expect the remaining excitation to propa- 
gate in its original spatiotemporal mode. For a coher- 
ent input pulse, one similarly expects the wavepacket to 
evolve with some probability into a correlated train of 
blockade-radius-separated photons. Furthermore, if the 
blockade radius is larger than the EIT-compressed pulse 
length, there may be a regime in which such a system 
could function as a deterministic single-photon source. 

In summary, we have shown that Rydberg blockade in 
EIT media can be harnessed for inducing strong photon- 
photon interactions, with applications to generating non- 
classical states of light, implementing nonlinear photonic 
gates, and studying many-body phenomena with strongly 
correlated light. This work opens several promising av- 
enues of research. With an eye towards single-photon 
generation, one can extend the presented wavefunction 
treatment to a density matrix approach and explicitly 
analyze the propagation of the remaining excitation af- 
ter the interaction- induced decay of multi-photon states. 
In addition, a gas of bosons (Rydberg polaritons) with 
a hard-sphere core (of radius Zb(B)/2) can be investi- 
gated both theoretically and experimentally in the co- 
propagating case. In particular, the previously neglected 
effects of es- endow these bosons with an effective mass 
oc —idj I \Lv g Y\ which plays a significant role for prop- 
agation distances larger than those considered in the 
present Letter. By including the effects of the coordi- 
nates transverse to the propagation axis, one can ex- 
tend this problem to higher dimensions. Furthermore, for 
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